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摘要 : 其 长 基线 干涉 测量 (very long baseline interferometry, VLBI) 是 重要 的 深 空 探测 器 测 轨 
HER. HE VLBI 网 (Chinese VLBI network, CVN) 在 中 国 探 月 工程 的 快速 精密 测定 轨 任 务 中 
发 挥 着 不 可 替代 的 作用 。 根 据 探 测 器 的 轨道 、 信 标 特征 和 数 传 信号 的 编码 、 调 制 方式 ， 研 究 数 传 
信和 号 的 仿真 方法 ， 并 仿真 上 升 器 月 面 起 飞 前 测 站 记录 的 基带 数据 。 仿 真 与 实测 数据 处 理 结果 的 对 
比 ， 表 明 该 仿真 方法 能 够 较为 真实 地 反映 数 传 信 号 对 DOR 侧 音 的 干扰 ， 且 能 够 对 测控 和 数 传 信 
号 同时 工作 的 VLBI 数据 处 理 系统 的 性 能 进行 检测 。 
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1 引 言 


其 长 基线 干涉 测量 (very long baseline interferometry, VLBI) 是 重要 的 深 空 探测 器 测 
轨 技 术 ， 具 有 空间 分 辨 率 高 、 测 量 精度 高 的 特点 。 中 国 VLBI 网 (Chinese VLBI network, 
CVN) 自 成 立 以 来 便 在 我 国 探 月 工程 的 快速 精密 测定 轨 任 务 中 发 挥 着 不 可 蔡 代 的 重要 作 
= j”"。 剥 娥 五 号 任务 是 我 国 探 月 工程 “ 绕 、 落 、 回 ”战略 的 第 三 步 ， 是 目前 中 国航 天 最 复 
© 杂 、 难 度 最 大 的 任务 之 一 ， 需 要 突破 月 面 采 样 、 月 面 起 飞 、 月 球 轨道 交会 对 接 和 地 球 大 气 高 
速 再 入 返回 等 关键 技术 。 探 测 器 由 轨道 器 、 返 回 器 、 着 陆 器 和 上 升 器 四 部 分 组 成 ， 在 整个 飞 
行 过 程 中 多 次 分 离 ， 存 在 四 器 组 合体 、 轨 返 组 合体 、 着 上 组 合体 和 单 器 等 工作 状态 ”。 然 而 
多 种 工作 状态 存在 X 频 段 多 频率 同时 工作 的 兼容 性 问题 。 如 在 上 升 器 月 面 起 飞 时 ， 着 陆 器 发 
送 的 数 传 信号 对 上 升 器 发 送 的 测控 信号 造成 干扰 甚至 淹没 部 分 频 点 。 
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四 是 月 面 起 飞 时 信 标 的 分 布 示意 图 。 上 升 器 有 两 组 测控 信 标 F3 和 F4。VLBI 观测 使 


的 差分 单程 测 距 (differential one-way ranging, DOR) 侧 音 为 点 频 信号 ， 通 过 调 相 的 方式 
调制 在 测控 信号 下 行 载波 上 。 在 XX 波段 上 ， 每 组 测控 信号 包含 两 对 DOR 侧 音 ， 以 F3 


EER 


为 例 ，F3PC 表示 下 行 信号 的 主 载波 ，F3L2, F3L1, F3R1 和 F3R2 依次 表示 与 主 载波 相距 


3.8, +3.8 和 


19.2, 


= 


+19.2 MHz 的 DOR. 侧 


Wr. 


音 ， 其 中 距离 主 载波 43.8 MHz 的 一 对 
位 测量 的 Qo 模糊 度 解 算 ， 距 离 主 载波 +19.2 MHz 的 一 对 用 于 VLBI 时 延 解 算 。 数 传 信 
使 用 二 元 相 移 键 控 (binary phase shift keying, BPSK) 调 


Hl Fz sk, 


ITAA 


其 频谱 扩展 范围 与 调制 码 


速率 成 正比 。 着 陆 器 的 数 传 信 标 为 F6， 码 速率 为 5 Mb/s， 其 频谱 扩展 范围 较 大 ， 会 对 邻近 
的 测控 信 标 造成 干扰 ， 进 而 使 VLBI 观测 结果 数据 质量 下 降 。 


VLBI 观测 结果 为 信号 同一 波 前 到 达 基 线 两 端 测 站 的 时 间 差 ， 简 称 时 延 ”。 
利用 时 延 和 测 距 测速 数据 实现 定 轨 定 位 和 轨道 


明确 


上 6 
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图 1 


HERT. E VLBI MEAS, MAAA 


上 升 器 月 面 起 飞 时 的 信 标 分 布 示 意图 


CVN 综合 


预报 ， 对 时 延 数据 的 实时 性 、 可 靠 性 和 精度 


机 对 测 站 记录 的 数据 做 互相 关 得 到 基 


线 上 的 互相 关 功 率 谱 ， 条 纹 拟 合 从 功率 谱 中 得 到 时 延 结 果 。 在 嫦娥 五 号 上 升 器 起 飞 前 ， 测 


= 所 .三 


控 信 号 受 


到 数 传 信号 的 干扰 ， 使 两 组 信 标 得 到 的 时 延 数据 质量 不 同 程度 地 下 降 。 


Pa 
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F6 的 包 络 遍布 在 观测 通道 内 ， 红 色 的 是 相关 相位 谱 ; 


横 轴 


是 观测 通道 内 的 频率 ， 通 道 带宽 


为 2 MHz， 频 谱 分 辨 紊 约 为 2 kHz. F3R2 距离 F6 的 中 心 频 率 约 12 MHz， 受 到 F6 第 二 旁 
六 的 干扰 ， 在 功率 谱 ( 见 图 2a)) 上 仍然 有 可 见 的 峰值 。 而 F4R2 距离 F6 的 中 心 频率 不 到 


3 MHz， 受 到 


余 时 延 


F4 的 结果 。 
受到 的 干扰 


在 
加 ， 然 


F6 主办 的 干扰 ， 其 功率 谱 ( 见 


， 横 轴 是 时 间 ， 纵 轴 


F3 受 数 传 信号 的 干扰 较 小 ， 


而 工作 频率 带宽 范围 


胶 大 ， 时 延 值 缺失 约 70%, eat 


图 2b)) 被 完全 淹没 。 
是 残余 时 延 值 ， 每 5 s 一 个 点 。 蓝 色 的 是 F3 的 结果 ， 红 色 的 是 
可 以 获得 不 间断 的 时 延 数 据 用 于 后 续 的 定 轨 ; F4 
ERIE 下降。 


TIR, VLBI XUY 


K B Æ BJ-TM 基线 的 残 


后 续 控 月 任务 中 ， 在 轨 探 测 器 数目 会 增多 ， 测 控 与 数 传 信号 同时 工作 的 时 间 段 会 增 
上 更 加 容易 出 现 多 频率 同时 工作 的 兼容 性 问题 。 
VLBI 信号 仿真 可 以 根据 目标 的 轨道 和 信 标 特征 生成 各 测 站 终端 记录 的 基带 数据 ， 常 
VLBI 数据 处 理 系 统 的 性 能 测试 。 仿 真 数 传 信号 可 月 


JF 


昌 于 数 传 和 测控 信号 同时 工作 时 的 系统 检 
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YE: a) F3R2 信 标 仍然 可 见 ; 


b) F4R2 信 标 被 淹没 。 
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3 ”上升 器 月 面 起 飞 前 后 BJ-TM 基线 的 残余 时 延 
测 ， 保 证 VLBI 信 标 在 实时 任务 中 的 可 靠 性 。 


ay VLBI 信和 号 仿真 软件 包括 Meyer-Zhao 等 人 开发 的 Datasim 和 Jakob 


的 VieRDS™ 


电源 信号 用 于 测 地 VLBI 项 目 。Datasim 首先 仿真 频 域 信 号 ， 多 


， 可 以 仿真 射电 


后 根据 观测 频率 通 道 选择 信号 变换 到 时 域 ， 忆 


再 根据 同一 波 前 


保存 为 VDIF 格式 。Datasim 可 生成 阿 塔 卡 玛 训 米 / 亚 训 米 波 阵列 望远镜 (Atacama Large 


Millimeter/submillimeter Array, ALMA) fi AXE 


率 的 基带 数据 ， 并 用 于 DiFX 的 性 能 


试 。VieRDS 全 部 在 时 域 上 操作 ， 


首先 仿真 基带 信号， 


然后 做 添加 时 延 模型 


! 等 操作 ， 


测 


输 


的 


VDIF 格式 数据 可 用 于 分 析 VLBI 全 球 观测 系统 (VLBI Global Observing System, VGOS) 


oe 


X 


宽带 观测 数据 处 理 链 中 的 系统 性 影响 。 
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本 文 根 据 探 测 器 的 轨道 、 信 标 特征 和 数 传 信号 的 编码 、 调 制 方式 ， 研 究 仿真 数 传 信 号 基 
带 数 据 的 方法 。 该 方法 根据 VLBI 系统 的 观测 流程 正 向 设计 ， 具 有 可 扩展 性 ， 满 足 多 目标 源 
音 号 仿真 的 需求 ， 且 全 部 在 时 域 上 操作 。 然 后 对 仿真 数据 解 调 、 解 码 ， 验 证 编码 、 调 制 方式 
实现 的 正确 性 。 最 后 本 文 根 据 嫦娥 五 号 的 实际 参数 ， 仿 真 月 面 起 飞 前 上 升 器 测控 信号 和 着 陆 
器 数 传 信号 同时 工作 的 VLBI 基带 数据 。 通 过 将 仿真 与 实测 数据 处 理 结果 进行 对 比 ， 表 明 该 
仿真 方法 能 够 较为 真实 地 反映 数 传 信号 对 DOR 侧 音 的 干扰 ， 且 能 够 在 任务 前 针对 测控 和 数 
传 信号 同时 工作 的 情况 检测 VLBI 数据 处 理 系 统 。 


2 数 传 信号 仿真 方法 与 验证 


在 深 空 通信 信道 中 ， 差 错 主要 有 两 类 ， 随 机 且 独 立 存在 的 随机 差错 和 密集 出 现 且 前 后 有 具 

了 相关 性 的 突 发 差错 。 为 保证 数据 传输 的 可 靠 性 ， 嫦 娥 五 号 探测 器 采用 国际 空间 数据 系统 咨 

询 委员 会 (Consultative Committee for Space Data Systems, CCSDS) 推荐 的 级 联 码 编码 和 

BPSK 调制 ”发 送信 息 数据 。 

Æ VLBI 系统 SLAB) 中 ， 探 测 器 发 出 的 射频 信号 被 多 个 测 站 接收 ， 并 与 氧 钟 提 供 的 本 

振 信 和 号 混 频 后 变 成 中 频 信号 ， 再 经 过 数据 采集 终端 滤波 、 下 变频 、 数 字 化 采样 和 量化 后 变 成 

基带 数据 ， 最 后 存储 在 磁盘 上 或 通过 高 速 数据 通信 网 络 实时 传送 到 VLBI 数据 处 理 中 心 。 在 

VLBI 数据 处 理 系统 中 ， 相 关 处 理 机 对 基带 数据 做 互相 关 处 理 ， 得 到 基线 上 的 干涉 条 纹 ; 后 
处 理 根据 干涉 条 纹 解 算 VLBI 时 延 。 


射频 信号 


磁盘 记录 
或 实时 传输 


E4 VLBI 系统 示意 图 


探测 器 信号 有 明确 的 时 域 表 达 形 式 。 从 探测 器 发 送 到 各 测 站 接收 所 经 历 的 光 行 时 能 
精确 反映 信号 同一 波 前 到 达 各 测 站 的 时 间 差 ， 因 此 仿真 数 传 信号 采用 的 时 延 模型 包括 光 行 
时 、 介 质 误差 和 仪器 误差 。 光 行 时 根据 探测 器 的 轨道 和 各 测 站 的 位 置 计算 ， 介 质 误 差 是 信号 
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在 大 气 和 电离 层 中 传播 产生 的 时 延 : 仪器 误差 包括 原子 钟 差 、 电 缆 时 延 和 观测 装置 内 部 时 
延 等 。 

仿真 射频 信号 的 数据 量 庞 大 。 嫦 娥 五 号 的 数 传 信号 在 XX 波段 ， 其 中 心 频率 约 为 
8.4 GHz。 根 据 奈 奎 斯 特 采样 定律 ， 仿 真 需要 16.8 GHz 以 上 的 采样 率 ， 难 以 工程 实现 。 射 
频 信和 号 与 本 振 信号 混 频 后 的 中 频 信号 通常 只 有 几 十 MHz 的 带宽 ， 且 包含 原本 的 信和 号 结构 及 
言 息 。 因 此 仿真 可 以 从 中 频 信号 开始 实现 。 

数 传 信号 仿真 根据 VLBI 系统 观测 流程 正 向 设计 ， 首 先 根 据 编 码 、 调 制 方式 和 时 延 模型 
构造 中 频 信号 ， 然 后 模拟 数据 采集 终端 滤波 、 下 变频 和 减 采 样 得 到 基带 信号 ， 最 后 量化 存储 
为 Mark5B 格式 。 第 世 I 节 介 绍 CCSDS 的 级 联 码 编码 方式 ， 第 区 3 节 讨 论 中 频数 传 信号 仿 
真 方法 ， 第 区 3 节 介绍 从 中 频 信 号 中 提取 基带 数据 的 方法 ， 第 区 引 节 介 绍 仿真 数 传 信号 的 验 
证 方法 。 
2.1 CCSDS 的 级 联 码 编码 方式 

CCSDS 推荐 以 里 德 - 索 罗 蒙 (Reed-Solomon, RS) 码 为 外 码 、 以 卷 积 码 为 内 码 的 级 联 
码 ”。 图 日 是 该 级 联 码 的 编码 框图 ， 其 中 RS 码 用 于 纠正 突 发 差错 ， 卷 积 码 用 于 纠正 随机 差 
并 ， 交 织 可 避免 在 连续 发 生 突 发 差错 时 超出 RS 码 的 纠 错 能 力 ， 伪 随机 化 使 用 收发 双方 约定 
的 伪 随 机 序列 加 密 通信 "。 


[RS 编码 -一 交织 上: 伪 随 机 化 H 卷 积 编码 ] 
图 5 ”级 联 码 编码 框图 


在 编码 代数 中 ， 信 息 序列 可 以 表示 成 多 项 式 ”: 
L-1 
m(x) = > mx , (1) 
i=0 


HP, L 为 序列 的 长 度 ，m; 为 第 i 个 码 元 的 值 。 

RS 码 的 码 元 是 伽 罗 华 域 (Galois field, GF) 上 的 元 素 ， 每 个 码 元 由 m bit 组 成 。GF(2™) 
上 的 全 部 元 素 由 一 个 m 阶 的 本 原 多 项 式 P(r) 推算 出 ， 包 括 27-1 个 非 0 元 素 {a} (i = 
0,1,2,.… ,2™-?) 和 一 个 0 元 素 ， 每 个 非 0 元 素 都 是 特殊 元 素 a Ne, a 是 本 原 多 项 式 的 
He. (n, k)RS 码 的 码 字 包括 n Em, HPE k MEREM n- k NARE. RS 编码 
首先 计算 生成 多 项 式 g(x): 


n—k 


g(t) = ][@- 4) . (2) 


i=1 


信息 多 项 式 m(z) 右 移 n 一 上 位 后 除 以 生成 多 项 式 g(z)， 其 余 式 为 监督 多 项 式 p(z): 


p(x) = x" mp(z)modg(z) , (3) 


其 中 ，mod 表示 求 余 式 。 因 此 对 信息 序列 做 RS 编码 的 数学 表达 式 为 ”: 


si(t) = p(z) +a" “m(z) . (4) 
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(255, 223)RS 码 的 每 个 码 元 由 8 bit 组 成 ， 其 本 原 多 项 式 为 F(x) = 28 +r +r? +241. 

交织 可 以 将 信道 中 产生 的 连续 突 发 差错 分 离 ， 避 免 超 出 RS 码 的 纠 错 能 力 ， 提 高 通信 系 
统 的 抗 干 扰 性 能 ”"。 交 织 是 将 原始 码 字 中 相 邻 的 码 元 分 离 ， 填 充 其 他 码 字 的 码 元 。 交 织 前 同 
一 码 字 中 相 邻 的 码 元 在 交织 后 的 距离 称 为 交织 深度 。 图 日 是 一 个 交织 深度 为 4 的 示例 ， 示 
例 中 有 4 个 原始 码 字 ABCD， 每 个 码 字 均 及 n 个 码 元 ， 以 A; 为 例 ， 它 表示 A 码 字 中 的 第 i 
个 码 元 。 通 过 交织 可 得 到 编码 信息 sz(z)。 


A, | B, | ©, 


D, : A, | B.| ¢,] D, 


6 ”交织 示例 (交织 深度 为 4) 


伪 随 机 化 使 用 收发 双方 约定 的 伪 随 机 序列 a(z)。CCSDS 推荐 的 生成 多 项 式 为 h(x) = 
z8 十 z7 十 z5 十 Z3 二 1” ， 其 阶 数 为 8。 序 列 前 8 位 初始 化 为 全 1 状态 ， 对 于 任意 aili > 8) A: 


8 
j=l 


当 生 成 多 项 式 h(x) 为 既 约 多 项 式 时 ， 伪 随机 序列 a(z) 具有 最 长 周期 ， 即 阶 数 为 8 时 ， 在 
255 位 后 重复 ”。 将 伪 随 机 序列 a(x) 与 交织 后 的 编码 信息 sz(z) 做 模 二 和 运算 ， 可 得 到 新 的 
编码 信息 53(7): 


$53(7) = a(x) © s2(z) . (6) 

卷 积 码 (K,1/n) 的 约束 长 度 为 KK ， 编 码 效率 为 1/n。 使 用 多 项 式 的 乘积 表示 时 域 卷 积 

运算 ，n 个 一 1 阶 生成 多 项 式 g(x) 构成 生成 矩阵 G(x) = [g1(7x),g2(x),… ,gn(x)]， 则 卷 
积 编码 可 表示 为 ”: 


sd4(Z) = 83(2)G(a) . (7) 
卷 积 码 (7, 1/2) 的 生成 多 项 式 为 gj(z) = 28 +15 +244 2941, ga(Z) = rê +r +r +r +l. 
2.2 ”构造 中 频数 传 信号 

数 传 信号 仿真 以 BPSK 信号 为 原始 信和 号: 


z(t) = s(t)V2P, sin(2rf.t) ， (8) 


RP, P, 为 探测 器 的 发 送 功率 ，f. 为 数 传 信号 的 调制 载波 ， 当 编码 信息 为 “0” 时 ，s(t) 取 
值 为 1， 对 应 初始 相位 为 0; 当 编 码 信息 为 “1” 时 ，s(t) 取 值 为 -1， 对 应 初始 相位 为 To 
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仿真 引入 时 延 模 型 d(t) 反映 数 传 信号 同一 波 前 到 达 不 同 测 站 的 时 间 差 。 由 于 探测 器 和 
各 测 站 的 相对 位 置 时 刻 变 化 ， 时 延 模 型 是 随时 间 变 化 的 连续 函数 ， 在 短 时 间 内 可 以 假设 为 线 
性 变化 ， 即 a(t) = T+ Tto WAE t 时 刻 接收 到 的 信号 实际 上 是 t 一 d(t) 时 刻 由 探测 器 发 
出 的 ， 因 此 测 站 接收 到 的 数 传 信号 为 : 


z(t) = s(t — To — nt) y 2P, sin [2r fe(1 — T )t — 2r feTo] ， (9) 
其 中 ， 书 为 测 站 接收 到 的 信号 功率 ， 计 算 方式 如 下 : 


P, = Ps + Po — Pa , (10) 


其 中 ，Pe 为 测 站 的 接收 天 线 增益 ，P 为 自由 空间 损耗 >。 对 比 式 (B) ARO aH RIK 
万 变 为 大 (一 五 )， 反 映 出 探测 器 和 测 站 的 相对 位 置 变化 导致 的 视 向 多 普 勒 频 移 。 

BPSK 信号 的 功率 谱 ( 见 图 团 在 频率 轴 上 以 调制 载波 大 为 中 心 向 两 边 无 限 延 伸 ， 其 主 
办 带宽 是 码 速 率 fo 的 两 倍 ， 旁 办 带宽 是 一 倍 。95% BPSK 信号 的 能 量 集 中 在 主办 和 第 一 和 
关上 。 为 降低 数据 处 理 所 需 的 采样 率 ， 仿 真 采 用 频率 fo = (1 — 71) fe — 2 fo 的 本 振 信号 与 测 
站 接收 到 的 射频 信号 做 混 频 ， 得 到 中 频数 传 信号 ， 表 示 如 下 : 


x3(t) = s(t — To — nt)v2P. sin (2n(fe(1 =n) — fo)t— Qn fero) ; (11) 


PE H E M PR 
图 7 BPSK 信号 功率 谱 
在 通信 链 路 中 的 诸多 噪声 也 是 仿真 的 重要 内 容 。 地 球 噪声 源 包 括 大 气 气体 的 辐射 、 云 二 
的 辐射 等 ， 地 球 外 噪声 源 包括 字 宙 背景 噪声 、 太 阳 和 月 球 辐射 、 天 体 射电 源 辐射 ， 人 为 噪声 
源 包 括 电子 电气 设备 的 无 意 辐射 、 输 电线 和 其 他 通信 系统 辐射 。 以 上 各 种 噪声 均 可 用 等 效 噪 
声 温度 衡量 。 噪 声 功率 P, 与 等 效 噪 声 温度 T, 的 关系 为 : 


P,=kBpT, . (12) 


其 中 , 大 为 玻 尔 效 曼 常数 ，Br 为 仿真 中 频 信号 带宽 。 仿 真 使 用 均值 为 零 的 高 斯 随机 序列 模 

拟 噪声 ， 其 方差 等 于 噪声 功率 。 中 频 仿真 信号 先 与 噪声 序列 辣 加 后 ， 再 提取 基带 数据 。 

2.3 ”提取 基带 数据 
VLBI 数据 采集 系统 对 中 频 信号 的 部 分 频率 滤波 、 下 变频 、 采 样 和 量化 存储 ， 为 VLBI 

数据 处 理 系统 提供 基带 数据 。 图 回 是 从 中 频 仿 真 信号 中 提取 基带 数据 的 流程 图 。 
仿真 首先 根据 探测 器 的 信 标 特征 设置 频率 通道 ， 并 根据 每 个 通道 的 天 空 频率 和 带宽 分 

别 设 计 带 通 滤波 器 。 滤 波 器 选择 凯 撤 窗 FIR 滤波 器 ， 它 具有 旁 久 抑制 性 能 好 且 相 频 特 性 呈 
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( Fes |} 通道 滤波 “上 和希 尔 伯 特 变换 | 人 ”下 变频 。 “] 
| 一 | ”量化 H 减 采样 一人 { 希 尔 伯 特 道 变换 ] 


8 ”提取 基带 数据 流程 图 


线性 的 特点 ”。 通 道 滤波 后 的 信号 是 实 信 号 ， 在 正 负 频 域 有 对 称 的 频率 分 量 。 为 避免 下 变频 
造成 正 负 频 域 的 信号 混合 ， 滤 波 后 的 信号 先 做 希 尔 伯 特 变换 只 保留 正 频 域 分 量 ， 青 做 下 交 
频 。 和 希 尔 伯 特 逆 变换 使 其 恢复 为 实 信号 。 中 频 信和 号 通过 减 采 样 降 为 基带 信号 。 量 化 采用 2 比 
特 量化 ， 输 出 二 进 制 结果 00, 01, 10 和 11 的 比例 分 别 为 18%, 32%, 32% 和 18%%， 此 时 量化 
效率 接近 最 大 值 " 。 量 化 后 的 基带 数据 以 Mark5B 格式 记录 在 磁盘 上 。 

2.4 ”仿真 数 传 信号 的 验证 

对 仿真 数 传 信号 解 调 、 解 码 ， 比 较 解码 信息 与 原始 信息 ， 可 验证 仿真 过 程 中 编码 和 调制 
方式 实现 的 正确 性 。 

BPSK 信号 解 调和 常用 相干 解 调 法 ， 其 关键 在 于 生成 一 个 同 频 同 相 的 相干 信号 。 中 频 仿 
真 数 传 信号 的 频率 为 fie = f-(l—71) 一 有， 相位 为 or = —27 fero M fon 表示 通道 的 天 空 
频率 ， 则 经 过 通道 滤波 和 下 变频 后 得 到 的 基带 信号 频率 为 fa = fe(l—71) 一 fm。 通道 滤波 
使 用 的 FIR 滤波 器 带 来 线性 的 相位 干扰 ， 即 中 频 信 号 在 提取 基带 数据 的 过 程 中 产生 一 个 固 
定 的 时 延 ， 该 时 延 可 表示 为 : 


Tar = 5 x AT ; (13) 


其 中 ，N 是 FIR 滤波 器 的 阶 数 ，AT 是 中 频 信号 的 采样 时 间 间隔 。 因 此 仿真 数 传 基带 信号 
解 调用 的 相干 信号 可 以 表示 为 ; 


c(t) = sin(2a fat + gm — 27 fiFTar) . (14) 


图 加 是 级 联 码 的 译 码 框图 ， 其 过 程 包 括 卷 积 码 译 码 、 解 伪 随 机 化 、 解 交织 和 RS 码 译 
码 。 卷 积 码 译 码 采 用 维特 比 译 码 算法 ”。RS 码 译 码 首先 求解 伴随 式 以 判断 是 否 存 在 误 码 。 
不 存在 误 码 时 ，RS 码 去 除 监 督 码 得 到 信息 码 ， 存 在 误 码 时 ， 译 码 通过 Berlekamp-Massey 
迭代 算法 和 Chien 搜索 法 得 到 错误 位 置 ， 再 通过 Forney 算法 计算 错误 图 样 ， 最 后 通过 错误 
图 样 纠正 信息 码 ”。 


卷 积 码 译 码 解 伪 随 机 化 


9 ”级 联 码 译 码 框图 


仿真 数 传 信号 经 过 相干 解 调和 级 联 码 译 码 后 得 到 解码 信息 。 在 没有 噪声 的 情况 下 ， 解 码 
信息 与 原始 信息 一 致 ， 表 明 仿真 数 传 信号 的 编码 、 调 制 方式 符合 预期 。 
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3 仿真 数 传 信号 对 DOR 信和 号 的 干扰 


2020 年 12 月 3 日 15 点 10 分 (UTC) MAS LA HARM MARK, BZ 
任务 中 最 重要 的 VLBI 观测 弧 段 之 一 。 然 而 在 起 飞 前 ， 着 陆 器 发 送 的 数 传 信号 对 上 升 器 的 
DOR 侧 音 造成 干扰 并 淹没 部 分 频 点 ， 使 VLBI 观测 数据 质量 下 降 。 图 中 展示 上 升 器 起 飞 时 
信 标 在 频 域 上 的 分 布 。 两 组 DOR 信和 号 与 数 传 信号 的 距离 不 同 ， 因 此 受到 不 同 程度 的 干扰 。 
F3R2 信 标 距离 F6 的 中 心 频 率 约 12 MHz， 在 互相 关 功 率 谱 中 仍 有 可 见 的 峰值 ， 条 纹 拟 合 
的 结果 虽然 存在 整 周 模糊 度 跳 变 ， 但 结果 完整 。F4R2 信 标 距离 F6 的 中 心 频 率 不 到 3 MHz, 
被 淹没 在 数 传 信 号 的 包 络 中 ， 没 有 条 纹 拟 合 的 结果 。 

VLBI 信和 号 仿真 参照 嫦娥 五 号 任务 中 的 预报 轨道 和 信 标 参数 ， 生 成 测控 与 数 传 信号 同时 
工作 时 测 站 记录 的 基带 数据 。 仿 真 观 测 时 间 从 当天 15 点 到 15 点 10 分 ， 为 月 面 起 飞 前 有 数 
传 信 号 干扰 的 时 段 。 观 测 信 标 在 X 波段 上 ， 其 中 上 升 器 发 送 两 组 测控 信 标 F3 和 F4， 其 等 
效 全 向 辐射 功率 (equivalent isotropically radiated power, EIRP) 为 -3 dB.W; 着 陆 器 发 送 
数 传 信 标 F6， 其 EIRP 为 25 dB.W， 码 速率 为 5 Mb/s。 测 站 有 BJ MTM, XH 16 通道 
记录 ， 其 通道 带宽 为 2 MHz， 使 用 2 bit 量化 。 仿 真 加 入 热 噪声 ， 其 中 宇宙 背景 噪声 2.7 K, 
月 球 噪声 290 K, BJ ATM 两 测 站 在 X 波段 上 的 系统 噪声 温度 分 别 为 60 K 和 30 K. 

仿真 数据 通过 相关 处 理 机 对 两 个 测 站 的 基带 数据 做 互相 关 ， 得 到 互相 关 功 率 谱 。 图 四 
是 仿真 信号 的 互相 关 功 率 谱 ， 蓝 色 的 是 幅度 谱 ， 红 色 的 是 相位 谱 ， 观 测 通道 带宽 为 2 MHz, 
谱 分 辨 率 约 为 2 KHz。 数 传 信 号 的 幅度 谱 以 调制 频率 为 中 心 向 两 边 延展 ， 其 包 包 络 遍布 在 测控 
言 号 的 观测 通道 内 。 在 观测 F3R2 的 通道 ( 见 图 10a)) P, DOR 信号 受到 的 干扰 较 小 ， 在 
F3R2 所 在 的 频 点 上 仍然 有 可 见 的 波峰 。 然 而 F4R2 的 功率 谱 ( 见 图 10b)) 被 淹没 在 数 传 信 
号 的 包 络 中 。 


20 50 
LF3R2 , 
ul alis i Addai ig 
= S A Wm 
= | Ses a Lie ef = Sa ay 3: ` 过 
HX A i ie ser El de Pe spi eae ages Fal 
E Bs a E T SA NA d HA 10 E g as cane } Pee n © 
| fe Peed, itd _10|。 到 ,各 Pid ii Se Ane 
l iad ae s N wS i 
tng 2 Pe 
—— -10 20 一 一 一 一 10 
= 0.5 1 0.5 2 00 0.5 1 0.5 2 
频率 /MHz 频率 /MHz 
a) b) 


YE: a) F3R2 信 标 仍然 可 见 ，b) F4R2 信 标 被 淹没 。 
10 “仿真 信号 BJ-TM 基线 的 互相 关 功 率 谱 


条 纹 拟 合 从 互相 关 功 率 谱 中 获得 VLBI NE. AOD 是 仿真 信号 的 残余 时 延 ， 每 5s 一 
个 值 ， 蓝 色 的 是 F3 的 结果 ， 红 色 的 是 F4 的 结果 。 由 于 仿真 使 用 的 探测 器 轨道 为 预报 轨道 ， 
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仿真 信号 的 残余 时 延 和 实测 数据 存在 一 些 差异 。 但 是 ， 对 于 被 淹没 部 分 频 点 的 F4 信 标 ， 数 


据 处 理 的 结果 都 表现 为 数据 量 大 幅 减少 ， 约 70%。 这 表明 数 传 信号 仿真 能 够 较为 真实 地 反 
映 数 传 信号 对 DOR 信号 的 干扰 情况 。 
e F3 © 
-2.5 e F4 
a my 
w] -3 a 2 > moe oS ee : 
= j EEE a a 
-3.5 i 4 = e 
15:00 15:02 15:04 15:06 15:08 15:10 
时 刻 
11 ”仿真 信号 BJ-TM 基线 的 残余 时 延 


在 探 月 工程 的 后 续 任务 
性 做 检查 。 若 发 现 问题 ， 
制 频 率 ， 使 数 传 信号 


信号 的 育 


则 调整 参数 予以 规避 ， 如 提高 


的 包 络 不 能 淹没 DOR 信 


前 可 以 使 用 设计 参数 仿真 数 传 信号 ， 对 VLBI 观测 信 标 的 兼容 
提高 DOR 信号 的 EIRP 值 


或 者 改变 数 传 


号 的 波 


标的 EIRP 后 的 仿真 数据 处 理 结果 ， 此 时 EIRP 为 15 dB.W。 从 互相 关 功 率 i 


We, KE 是 通过 提高 F4 信 


普 ( 见 图 12a)) 


中 可 以 看 到 DOR 信号 的 波峰 没有 被 数 传 信号 的 包 络 淹没 ， 因 此 条 纹 拟 合 得 到 的 残余 时 延 
( 见 图 12b)) 没有 大 量 缺 失 。 仿 真 数 传 信号 不 仅 可 以 在 后 续 任 务 前 规避 一 些 信 标 兼容 性 问题 ， 


而 且 可 用 于 辅助 开发 测控 信号 受到 干扰 时 的 数据 处 理 新 算法 。 
25 š 40 -2 
1 F4R2 [ * F4] 
-2.5 E 
fanl a .435 = n 
= 7 2 £ 
Hin Be TES | N a 8 3 
w 5p iy cel eal 到 E ee rs 
4 = oon wae EE nes oe {ja -3.5 Lon eT 
-5 ave | 
0 0.5 1 o 525 475.02 15:06 15:10 
频率 /MHz 时 刻 
a) b) 


VE: a) 互相 关 功 率 谱 ，b) 残余 时 延 。 


图 12 


提高 DOR 发 射 功率 的 仿真 数据 处 理 结果 
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4 结 论 


本 文 根 据 VLBI 观测 和 数据 记录 流程 研究 探测 器 数 传 信号 的 仿真 方法 。 该 方法 包括 信 
县 编码 调制 、 添 加 时 延 模型 、 构 造 中 频 信号 和 提取 基带 数据 等 步骤 ， 根 据 探测 器 的 轨道 、 信 
标 特征 和 数 传 信号 的 编码 、 调 制 方式 生成 各 测 站 记录 基带 数据 。 在 后 续 探 月 任务 前 ，VLBI 
信号 仿真 检测 数据 处 理 系统 不 能 忽视 数 传 信号 的 干扰 。 本 文 模 拟 嫦娥 五 号 任务 中 月 面 起 飞 
时 的 观测 场景 ， 生 成 BJ 和 TM 测 站 记录 的 基带 数据 。 仿 真 信号 的 互相 关 功 率 谱 表 明 ， 数 传 
信号 对 测控 信号 各 侧 音 造成 不 同 程度 的 干扰 ， 其 中 F3R2 信 标 部 分 受到 干扰 而 F4R2 信 标 被 
淹没 ， 进 而 F4 信 标 的 后 处 理 时 延 值 数据 量 比 F3 大 幅 减 少 。 该 结果 与 实测 数据 处 理 结果 具 
一 致 性 ， 表 明 该 仿真 方法 能 够 较为 真实 地 反映 数 传 信号 对 DOR 侧 音 的 干扰 ， 且 能 够 用 于 
在 后 续 探测 任务 中 的 VLBI 数据 处 理 系 统 的 性 能 检验 。 
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Abstract: Very Long Baseline Interferometry (VLBI) is an important technique for deep 
space probe orbit measurement. Chinese VLBI Network (CVN) plays an indispensable role 
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in the fast orbit measurement and determination of the Chinese Lunar Exploration Program. 
In the Chang’E-5 (CE5) mission, VLBI orbit measurement system utilizes Differential One- 
way Ranging (DOR) sidetone modulated on the probe downlink TT&C signal as beacon. 
Before the ascender takes off from the lunar surface, the DOR sidetone is interfered by the 
data transmission. Some frequencies are submerged. As a result, the observation quality of 
the corresponding beacons are greatly reduced. In the follow-up lunar exploration mission, 
the number of on orbit probes and the duration in which TT&C and data transmission signal 
work simultaneously will increase. However, the available frequency range is limited. The 
problem of interference between multiple probes and frequencies becomes more serious. In 
this work, we carry out study on the method of simulating data transmission signal, generate 
the baseband data recorded by the station before the take-off from the lunar surface of CE5, 
and compare the data processing results of simulation and actual observation. Our results 
demonstrate that the simulation truly reflects the interference of the data transmission signal 
to the DOR sidetone, which makes it possible to carry out test for the scenario of which 


TT&C and data transmission signal work at the same time. 
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